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This paper presents a general averaging procedure for a set of observers which are tilted with 
respect to the cosmological matter fluid. After giving the full set of equations describing the local 
dynamics, we define the averaging procedure and apply it to the scalar parts of Einstein's field 
equations. In addition to the standard backreaction, new terms appear that account for the effect 
of the peculiar velocity of the matter fluid as well as the possible effect of a shift in the coordinate 
system. 
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The standard cosmological paradigm implicitly relies on a spatial averaging procedure applied to both the space- 
' time metric and the Einstein field equations. This procedure is required to obtain an effective homogeneous model for 
, the large scale properties of the Universe, given a clumpy local distribution of matter that is known to emerge from 
P_, ■ gravitational accretion. In the context of perturbation theory, this effective homogeneous model is referred to as the 
O background space-time. Nevertheless, the late time Universe appears highly structured on scales up to a few hundreds 
of megaparsecs, but a scale of homogeneity seems to exist at larger scales p^]. Although this homogeneity scale arises 
$_i ' from observations made along the past lightcone of the observer, it is usually used to support the cosmological 
principle, that states the existence of a homogeneity scale for the spatial distribution of matter. It is beyond the scope 
of this paper to address the problem of averaging over the lightcone; we will rather stick to the conservative point of 
^SJ ■ view and suppose that there exists a spatial scale of homogeneity. The usual way to deal with this issue is to ignore 
J> ' the structures on all scales and consider the spatially locally homogeneous and isotropic solutions of Einstein field 
0^ , equations known as Friedmann-Lemaitre-Robertson- Walker (FLRW hereafter) metrics. To account for the presence of 
structures, the model then relies on a perturbative expansion of the metric components around this FLRW background. 
Even though this model has proved to be extraordinarily successful in explaining a lot of features of the observed 
Universe, through the so-called concordance model, it suffers from a lack of theoretical basis from the point of view of 
the averaging it involves. Indeed, one would like to see the background emerge from the averaging procedure, rather 
than being postulated from the very beginning. Moreover, it is known that the nonlinearity of Einstein field equations 
implies that the equations for the smoothed metric are not equivalent to the smoothed equations for the actual local 
metric of space-time 0. In other words, the coarse-graining of the local fluctuations influences the kinematics of the 
' effective model on large scales. This has been studied in the last decade for a simple averaging procedure, leading 
k>( ' to the concern with the backreaction effect [1, 01 . From a more observational point of view, this topic has recently 
' received a revival of interest as a possible natural explanation for the Dark Energy problem (see e^., for good reviews, 
I d, @] , and for recent attempts to obtain observational constraints on the homogeneous model 0, U ) that tarnishes the 
success of the concordance model (in the sense that this latter fails to explain elegantly the origin of the phenomenon). 
In this paper, we generalize the work done by T. Buchert about the backreaction effect for a set of observers comoving 
with the matter fluid by extending the averaging procedure to an arbitrary set of observers, hoping in particular 
that it could shed light on the gauge issue in the backreaction context (i.e. the gauge dependence of the backreaction 
effect in perturbation theory). Buchert formalism has been used successfully in the synchronous gauge [9], where it 
applies directly without ambiguity. But in the Newtonian gauge, for example, different authors disagree on the type 
of effect obtained [a[iil,[il. These differences can be traced back to the fact that the different authors have used 
different definitions of the effective Hubble parameter. In our work, we define it with respect to the expansion of 
the cosmic fluid as seen by the observers. It seems that this latter quantity is more suitable to describe a physical 
expansion as well as to compare the gauges between each other. We aim at studying in detail the Newtonian gauge 
in this context in a forthcoming paper [l^. Moreover, a recent model developed by Wiltshire suggests that the 
choice of the set of observers relative to which the averaging is performed is a crucial feature. 
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The paper is organized as follows. Section II presents the formalism used to address the problem; we insist in particular 
on the foliation of space-time by a set of observers that are not comoving with the matter flow, and give the general 
equations governing the dynamics of both the geometry and the quantities associated with matter. In Sec. Ill, we will 
present the generalized averaging procedure that aims at keeping the main characteristics of the one introduced by 
T. Buchert. This will require a discussion about what really has to be average in order to define an effective Hubble 
parameter for the model. Thanks to this way of averaging, we will then give explicitly the set of equations that 
governs the kinematics of the averaged model, and comment on the new terms arising from the tilt between observers 
and the matter flow. Finally, we sum up the results and discuss possible and ongoing uses that can be made of this 
new formalism in Sec. 4. 



II. SPACE-TIME DYNAMICS IN AN ARBITRARY COORDINATE SYSTEM 

The cosmological model consists in a space-time M endowed with a metric tensor gab, and filled with a matter 
content characterized by its energy-momentum tensor. We will assume that gravitation is well described on all scales 
by general relativity, so that the metric and the matter content are linked by the Einstein field equations. The line 
element can be written: 

ds^ = gabclx^dx^ , (1) 

where (a;")ae{o,i,2.3} are the space-time coordinates. Throughout the paper, we will employ Latin letters of the 
beginning of the alphabet (a, b,c, ...,h) to denote space-time indices, and Latin letters of the end of the alphabet to 
denote spatial indices (this concept will be made clearer in the next section after the projector on spatial hypersurfaces 
has been introduced). 

A. Foliation of space-time 

Let us start by setting the notations employed to describe the 1-1-3 foliation of space-time that will be used in this 
paper. It consists of introducing a set of observers 0{p) defined at each point of the space-time manifold p ^ A4, 
characterized by a unit four- velocity field n° that is everywhere timelike and future directed, i.e. n^ria — —1. Once 
this four-velocity is known at each point of the manifold, one can foliate space-time by a continuous series of spacelike 
hypersurfaces that are simply the hypersurfaces everywhere orthogonal to the vector field n". Moreover, after having 
defined the projection tensor field on these hypersurfaces by hab '■= gab + naUb, one requires the four-velocity of 
observers to be irrotational, i.e. h'^h'^n^a-h^ = 0; this last property ensures that the projection tensor hab is also a 
well-defined Riemannian metric tensor for the hypersurfaces orthogonal to the four- velocity n". The line element can 
then be written, with respect to this foliation'^: 

ds^ = -{N^ - N,N^)dt^ + 2N,dtdx' + h.^dx'dx^ , (2) 

where we have introduced respectively the lapse function N{x°') and the shift three- vector N^{x"') such that the 
four- velocity of observers can be decomposed as follows: 

n'' = i^(l,-7V*) ,na = iV(-l, 0,0,0) . (3) 

Two other quantities characterizing the hypersurfaces orthogonal to n° will be very important in the following: 

• the intrinsic curvature of the hypersurfaces: TZ := h'^^TZab, where TZab is the three- Ricci curvature of the 
hypersurfaces, 

• the extrinsic curvature (or second fundamental form): Kab ■— ~h'ahfnc-,i that encodes the way the hypersurfaces 
are embedded in the manifold M.. 



^ For a clear description of the 1+3 foliation of space-time, see, e.g., the very clear paper by Smarr & York (in particular for a physical 
interpretation of the lapse and shift) [Tsll and more recently, the lecture notes [161 . 
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Up to now, we have introduced a set of observers and the purely geometrical quantities associated with it. In order 
to have a complete description of the cosmological space-time, one still has to prescribe the matter content of the 
Universe. In the following, we will assume that this matter content can be well described by a perfect fluid (not 
necessarily irrotational) of energy density pressure p{x°') and 4-velocity u"-{x^) (with Uau"- = —1), so that its 

stress-energy tensor reads: 

Tab = (P + P)UaUb + PQab ■ (4) 

Note that in this work, the four-velocity of matter m° is not necessarily aligned with the four-velocity of the observers 
n", so that one can introduce a vector field w° corresponding to the relative velocity of the matter fluid with respect 
to the fundamental observers, is spacelike and orthogonal to n° {v°-na — 0) and one has: 

= 7(n'^ + V-) with 7 = ^ , (5) 

where 7 is the usual Lorentz factor (see in particular [l7j for a 1-1-3 covariant treatment of multifluids cosmology, in 
which the effect of the tilt between a fluid and the observers is described). 

We now have all the ingredients to write down the Einstein field equations. Thanks to the 1+3 foliation, they can 
be separated in two different sets: the constraint equations that have to be satisfied on every hypersurface, and the 
evolution equations, that prescribe how the fields {hab,Kab) evolve from one hypersurface to another infinitesimally 
closed. The so-called Hamiltonian constraint reads: 

7^ - K}Ki = WnGe + 2A , e ^ Tabn^n" = 7^ + (7^ " 1)P , (6) 

where K := Kl. The momentum constraint is: 

V^K] - V,K - SvrGJ, , J, = -Tabu'^h] = j^ip + p)vj , (7) 

where we have defined the projected covariant three-derivative on the spatial hypersurfaces of any tensor field t"'" '^ ^ f- 
Vfji^ ' ^^a g := h'^h^, ...h'^ih'^ Vet"" '"'^a" c"- The evolution equation for the first fundamental form reads: 

^dth, = + ^V(,iV,) , (8) 

and for the second fundamental form, one has: 



^dtK] = U) + KK] - as; - ^V.V'TV + 1 (^KlVjN" - K^^^kN' -f N^\IkK'^ - SttG + ^(e - 3^)6]^ , (9) 



with Sij = 'j'^pViVj +p{hij +j'^ViVj). These equations have to be supplemented by the energy-momentum conservation 
for the matter fluid: aT^^ — 0. We will now introduce the standard decomposition for the covariant spatial derivatives 
of the four-vectors in terms of their trace, symmetric trace-free and antisymmetric parts. Writing / = n"Va/ for any 
quantity /, one can write: 

VaTlb = -UaUb + ^^/iab + ^ab (10) 

with ^ = hlh'^"'\/ and = hlh^S/^^nd) - ^C^afc ; 

VaUfc = -Ha {2hb + VVUb) ~ j\7 aVUb + ^9hab + (7ab + ^^ab (H) 

with 9 = hlh'^°-\IcUd , (^ab = h^h'^S/^cUd} - ^Ohab and uJab = hlhf,\I[cUd] ; 

VaWfc = -na{vb+ Xb) +Yanb+'^nhab + Pab + Wab (12) 



K = hlh''-°-\/ cVd , f3ab = h^^hfS/ (^cVd) - ^nhab 



with 

and Xa = n^KVcVb , n = n^hlVcVa , Wab = KhtV^^Vd] 



In these relations, ^, and k denote the isotropic expansion rates of the four-velocities n", u°' and of the peculiar 
velocity d° respectively, and Safe, ^^ab and (3ab their shears, with respect to the threading of space-time induced by the 
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vector n°. ujab and Wab are, respectively, the vorticities^ of and v°- in this same fohation. These quantities are those 
measured by the observers with four-velocities in their instantaneous rest-frame. In particular 0, aab and tOab differ 
from the usual expansion, shear and vorticity of the matter fluid as measured by observers comoving with this matter 
fluid (by acceleration terms essentially), that are defined by the decomposition of -t- u°-u'^){g^'^ + u^u'^V cUd)- For 
example, the expansions are linked by the relation: 

e = Va^i" = ^^ + 7(72^;%, . (13) 

Using ([5]), one can relate this quantities as follows: 

i = 7-16*- (14) 

Sq6 = l^^Oab - Pab - 7^ {^{ab) - ^Bhal}j (15) 
Wab = I'^^ab - l^B[ab] , (16) 

where we have introduced the tensor: 

Bab = ^l^iVarib + VaVb) + PcaV^Tlb + PcaV^Vb + WcaV^Hb + WcaV^Vb , (17) 

whose trace is given by i? = ^'^'^^ + PabV^'v^ , which encompasses the peculiar velocity effects. In our notation, angular 
and round brackets denote the antisymmetric and symmetric parts, respectively, of a tensor projected with hab- It is 
important to note that the extrinsic curvature Kij is related to the derivative of by: Kij = —hfh'j'Vanb, so that 
the equations for the expansion and shear can be obtained by suitably replacing Kij in the Einstein field equations 
(O-®. Let us finally introduce the following notation for convenience: 

Ob = --/K-j^B (18) 

<JBzj = -iPij -1^ {B(ij) -^Bh^j] , (19) 



so that: 



^ = ^-'{9 + 0b), (20) 

= 7"^(o'ij + o-Bi]) ■ (21) 

B. Field equations 

With these notations in hand, we can now write down the complete system of equations describing the behavior of 
matter on local scales. They read: 

-f^n + ^6*2 + + Uob - 2(7^ - 2(t| = {p + v^p) + 2A72 (22) 

^{\/,0 + ^jOb) ~ ^^a] - V,a^^- - « 0(0 + 0b)V,« + {o] + ct*^)V,«^ = 8^G7^(p + p)vj (23) 



^ Since it defines the foliation, is vorticity free by definition. 
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for the constraint equations ©-((Tl). Moreover, the evolution equations deduced from (l8])-(l9]) are: 

J^dthj + OB)h^J + 2((jy + asy) + ^V(jA^,) (24) 

7 

+7V^ViiV + N'^Vkid + Ob) + 47rG (iV7^p(3 - v'^) - N"fp{3 + 2j^v^)) (25) 
3 



^t^l = ^^—^ — -^^) + '^B,) - diOB, - Ni{n) - -7^<5}) - -{e + eB){cT] + a),^) 
+7(v*v,7v - iv"Vfe(7v)5;-) - VfcTvx^' + 4,0 



+V,N\al + a'sk) + A^'V^a} + iV^V^a^, + 8^G7V7'(p + - y^}) , (26) 

with (7^ = ^j'^i ^'^"^ '^s ^ '^bj'^bi/'^'^'^ij'^B ' FinaUy,the equations obtained from the energy-momentum conservation 
VaT° = read: 

-9tp+(«'=-— )Vfcp+(7-^e + 2i;« + «fca'=)(p+p) + -9tp + i;2 ^'^-—j V.p = (27) 

l5,.'=-^V..'=+(i^0-..a^)z;'= + a^+7-M4'-^+'-M + :^^ = 0. (28) 

We have noted Ofc := hak'h"' = {'^kN)/N the acceleration of the observers along their world-lines. 

III. THE GENERAL SPATIAL AVERAGING PROCEDURE 

In this paper, we will retain as our guideline the averaging formalism introduced in cosmology by T. Buchert 0, 
Namely, we will consider the operator acting on scalar quantities '>p{t,x^) such that, for any compact domain V 
included in an hypersurface orthogonal to the observers worldlines: 



(^)^= ^ / i;{t,x')Jd'x , (29) 

IV 



Vv 



where J = y^det {hi j) and Vt> = J Jd?x is the Riemannian volume of V. From Eq. ([24]) . one can deduce an evolution 
equation for the quantity J: 

IjdtJ ^-i~^N{9- k)+VuN'' , (30) 
that translates into a evolution equation for the volume Vp: 

dtVv = j {l^'^N{e - k) + VkN^^ JcPx . (31) 

One can immediately see that the volume scale factor as defined in the usual Buchert framework, by = 

1/3 

{VT>{t)/Vxi{ti)) is no longer related in a simple way to the averaged expansion of the matter flow as measured 
by the observers, {N9)i~,. But the measurable quantity of interest in cosmology is this latter expansion, not the 
scale factor, so in building an effective homogeneous cosmological model, we should concentrate on estimating the 
observable quantity, i.e. the expansion of matter. In the following, we will then retain as a definition for the Hubble 
rate of expansion: 



1 

Vv 



iH-D = {N9)^ = — / NO J (fix , (32) 



•D 



and define the effective scale factor ax>(t) such that: 



i/^ = ^ . (33) 
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The relation between this effective (matter expansion) scale factor and the (observers) volume scale factor is then: 

a^^apexp^^ (^(^{-/-^ ~ l)Ne ~ Nk + V kN'''^ {t )dt'^ . (34) 

The definition ((33|) coincides with the volume scale factor for observers comoving with the fluid (i.e. implying 7 = 1 
and K = 0) in a coordinate system comoving with the fluid (i.e. iV* = 0). This scale factor has been introduced 
only to rewrite the averaged equations in a form that can be compared with the usual Friedmann equations, but it 
is important to insist on the fact that everything could be expressed in terms of the effective Hubble expansion only, 
without any reference to a scale factor of an y k ind. It has to be noted that the observers volume scale factor a^J is 
the one whose evolution has been studied in |12| : it seems more natural to concentrate on the matter expansion rate 
as seen by observers, rather than to the expansion rate of the observers themselves. One of the main relations in 
averaged cosmologies comes from the fact that averaging and evolution do not commute, which, making use of Eq. 
([?T|) translates into: 

[dt., (.)i,]V'(i, ^1 = {{N7-\e + 0b) + VfcA^^-) - (^N^-\e + Ob) + ^kN'')^ , (35) 

for any scalar function tp{t, x^). One can then write the averaged equations arising from averaging the scalar parts of 
Einstein field equations, Eqs. (HH) and (PSj) : 

6Hl = IGttG {{l^N^p)^ + {iHl^ - + 2A (iV^')^ - (7'^^'7^)^ - Qv + Cv (36) 



•^'"^ = -47rG(iV^(7V(l + ^'') + 3p(l + 27V)))^+A(7V2^)^ 



+Qv + Vv+JCv+Tv- Cv , (37) 
where we have defined the standard kinematical backreaction: 

Qv^l {{{N9f)v - ml) - 2 {N^^^)r, . (38) 
and additional backreaction terms as: 

Cv ^ 2 {N^al)^ - \ {{NObY)^ - \ {N^OOb)^ (39) 



Vv = (^atiV)p + <j7iVVfcV^A^)^ (40) 



I /„, AfV7 . Yvfc Aj\ 

JCv = {NN''\7kO)^ + (^NN^^/kOB)^ + ((A^7"'^?b + VkN^)Ne'^^ - 3 (^N-f-^eg + ^kN^}^ Hv 

- {NdteB)v + {N^l^^e)^ + {N'j-^Ob)^ - 2 (N^el)^ (41) 

-I {N'eli^-' !))„ - I {N^eOBij-' - !))„ - 2 (N^alil-' l))v • (42) 

Note that for the sake of simplicity, we have used 7 = n°-V a'^ — {dt^)/N — {N'^V i^) / N . 

Averaging the energy conservation equation (|27p in the same way, one can also obtain the equation for a scaled 
averaged energy density N'^ p: 



dt (N^p)^ + 3Hv {{N^p)^ + {N'p)^) - 2 (pNdtN)^ - [9*., {.)v]iN'p) + [{Nj-'dB + VkN'')N\p + p)'^ 

- {n^-'9b + VkN'^)^ {{N'p)^ + (N'p)^) - {N9ij-' 1))^ {(N^p)^ + (N^p)^) (43) 

- l^N^v'' ~ iV'=)Vfe(p + p))^ - {N\2vv + Vka''){p + p))^ - (N^^dtp)^ . 
Finally, one can obtain an integrability condition for the system formed by Eqs. p6p and (|37p : 

dtQv - dtCv + dt (7'iV'7^)^ - 2Aat (N^j^)^ + 2Hv {Qv - SCv + {I'^N^U)^ + 2Vv + 2ICv + 2Tv) 

= IGttG {dt {i^N^p)t, + Hv [{N^P-/^ (27 + 1 + v^))^ + {N^p {j^v^ + 37(1 + 2-/^v^)))^] + dt {j^v^ N^p)^]{U) 
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Unfortunately, when v'^ is not zero, this last equation does not simplify because p and p appear with different space- 
time dependent coefficients in ([55]) and (|37p . so that one cannot make use of P5|) . Nevertheless, it is easy to see that 
when = and N*" = 0, the system ([36])-(l37])-(l43l)-(|44|) reduces to the well known results of i-e. that the terms 
Ct>, ICt> and J-p identically vanish. Another remark on this system: one can note that Qt> and Ct> appear everywhere 
in this system in the combination Q-p — Ct>, so that in a tilted coordinate system, it is more natural to define the 
kinematical backreaction as Qv — .Cv, Cv taking into account the pure velocity effects; nevertheless, we chose to stick 
to the conventions of [3] so that the contact with the comoving observers can be done easily. First, if applied to the 
Newtonian gauge of perturbation theory, this formalism is expected to differ from previous analyses, in particular the 
one of [Hj, because of the additional terms L-d, A^d and JFp. Again, this comes from the fact that they defined the 
Hubble parameter with the expansion of the observers rather than with the expansion of the matter fluid. Second, 
the model proposed in [3] relies mainly on the fact that the lapse function is not constant throughout space-time due 
to the status of observers in virialized regions with respect to averaged cosmological observers identified to be void 
observers. But it is clear that if the void observers define the normal vector n", then the virialized observers must 
have a peculiar velocity with respect to them, so that the new terms also have to be estimated in this context if one 
wants to have a fully consistent framework. 

IV. CONCLUSION AND OUTLOOK 

In this paper, we have generalized the averaging formalism of [1, Q to account for a possible tilt between the 
cosmological matter flow and the set of observers for which the (spatial) cosmological principle is supposed to hold. 
We have shown that, even if the general picture is not dramatically changed, the averaged system of equations 
governing the behavior of the effective model on large scales involves a lot of nontrivial new terms that describe the 
velocity effects (as well as the effect of the presence of a shift vector in the metric). Our formalism, based on the 
standard Riemannian averaging, is coordinate dependent, as pointed out in a recent paper [isj . In this paper, the 
authors have shown that, in perturbation theory, the crucial point about the gauge dependence of the formalism lies in 
the change in the domain of averaging and they have explicitly built an alternative averaging procedure that is gauge 
independent. Therefore, in order to understand the physical meaning of the gauge dependence of our formalism, it will 
be employed in forthcoming works to study the backreaction from perturbations in the Newtonian gauge, and then 
to analyze the differences in the backreaction effect when one try and evaluate it in different perturbative gauges. 
A comparison of the results in the different gauges with the ones of the gauge independent procedure of [l8| also 
seems important. Moreover, the formalism presented here will be at the heart of an investigation of the link between 
lightcone and spatial averagings, a link that is necessary to put observational constraints on averaged cosmological 
models. Finally, the author thinks that this will be useful in describing more precisely Wiltshire's idea about the 
difference between cosmological observers and matter observers in virialized regions 
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